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Abstract 

Classical Gaussian white noise in communications and signal processing is viewed as the limit 
of zero mean second order Gaussian processes with a compactly supported flat spectral density as 
the support goes to infinity. The difficulty of developing a theory to deal with nonlinear transfor- 
mations of white noise has been to interpret the corresponding limits. In this paper we show that 
a renormalization and centering of powers of band-limited Gaussian processes is Gaussian white 
noise and as a consequence, homogeneous polynomials under suitable renormalization remain white 
noises. 
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1 Introduction 



White noise is critical in the development of statistical signal processing and in models for commu- 
nication channels. The centrality of the process is that any second order covariance function can be 
realized as the result of white noise through a linear filter. In the stationary case white noise is the 
basic building block of of constructing optimal filters. In the classical context, white noise is viewed as 
a limit of a second order process that has a flat spectral density of compact or finite support (referred 
to as the bandwidth of the process) as the support becomes infinite [221 HI E]. Such processes cannot 
be physically realizable because they would have infinite energy, and yet white noise plays a crucial 
role in developing practical filters. 

The difficulty of defining white noise is not just because of the infinite energy. Indeed defining a 
continuous-time white noise process presents difficulty even in the probabilistic context. This is because 
the induced probability measure cannot be countably additive on the space ^[0, T] i.e. constructing 
a Gaussian process whose covariance function is a Dirac delta function results in the underlying 
probability measure being only finitely additive, and thus the white noise map ntioj) is not a bona fide 
random variable [6| 1201 [TO] . Indeed because of the difficulty of mathematically dealing with white noise, 
Balakrishnan in a series of papers [U (5J [7] developed a finitely additive framework for analyzing white 
noise processes and the associated calculus. This was further developed through the idea of liftings 
in the work of Kallianpur and Karandikar |17j . An alternative approach exploiting the structure of 
abstract Wiener spaces can be found in the work of Kuo [18], and Gross [131 E] f° r example, where 
the idea is to work through on the lifted space through the lifting map that results working on the 
abstract Wiener space. 

Classical white noise is only defined as a generalized process, i.e., it does not induce a countably 
additive measure on the Hilbert space. However certain transformations of white noise do induce 
countable additive measures. This is closely related to the result of Sazonov [20] on the existence of 
countable additive measures on Hilbert spaces. In the linear context, any continuous linear operator 
such as a kernal operator acting on white noise results in the resulting map defining a bona fide 
stochastic process. The integral operator Ln(t) = f^n s ds as a mapping from L2[0,T] — > L2[0,T] 
induces a countably additive measure whose extension to C(0, T) [8J is the Wiener measure [TBI 
121] , Viewed this way, Gaussian white noise is the derivative of the Wiener process even though 
the Wiener process is not differentiable almost surely for any t. However when restricting ourselves 
to linear operators we can make sense because in most applications (filtering and communications) 
such operators are Hilbert- Schmidt [18] and we can interpret the results both probabilistically as 
well as algorithmically. For nonlinear transformations the situation is more complicated. Indeed one 
class of nonlinear transformations that induce countably additive measures is that associated with 
S-continuous mappings where the S-topology is that which is associated with semi-norms associated 
with nuclear operators [21 El [9] . 

There still remains the question whether nonlinear transformations, the powers for example, of 
white noise can be suitably interpreted, even as generalized processes 110] . The need for such 
interpretation can be found in many applications in mathematical and quantum physics [U [19] for 
example. Indeed, in [19] the authors provide a heuristic justification for treating a renormalization of 
a squared white noise term as white noise. 

In this paper we show that under a suitable renormalization, integral powers of Gaussian white 
noise viewed as the limit of a band-limited Gaussian process with flat spectral density is indeed 
Gaussian white noise, a non-trivial fact given that non-linear transformations of Gaussian random 
variables are not Gaussian. 
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2 Some preliminaries 

Let us first see some of the classical results related to white noise. For simplicity we restrict ourselves 
to real-valued processes. The extension to 3ft d valued processes is direct. We denote the inner product 
in L 2 [0,T] by [f,g] T = f f{s)g(s)ds, f,g G L 2 [0,T] and the norm is denoted by [f,f] T = WfWx- The 
space L 2 (— 00,00) is simply denoted as L 2 

Let {X\y(t), —00 < t < 00} denote a zero mean stationary Gaussian process. Let R(t) = 
JE[X(t + s)X(s) denote the covariance and it is assumed that \R(t)\dt < 00. By Bochner's 
theorem there exists a spectral density 5(A), —00 < A, 00 and 5(A) = R(t)e~' l27rXt dt. 

Now suppose 5(A) is band- limited and flat as follows: 

S w (\) = 1 , —W < X <W 
= otherwise 

Now from the fact that 5(A) has support in [-W, W] it follows that there exists a spectral process 
that is Gaussian and independent on non-overlapping intervals denoted by Xw(\) such that for every 
(a, b), j"s w (X)dX = nJaJaX(d\)X(dX')\ and 

/oo 

e i2nXt X w (X)dX, -00 < t < 00 (2.1) 
-00 

Moreover the limit in mean square (denoted by q.m) 

lim X w (t) q ' m X(t) (2.2) 

W— >oo = 

exists and is called Gaussian white noise. 

In particular we see that such a limiting process process will have the following properties: 

i) The random variables y = Jq X(s)f(s)ds will be distributed as N(0, where / G L2[0,T]. 

n)TE([f,X] T [g,X] T ) = [f,g] T 

hi) Let (f>i(t) be an orthonormal functions in L 2 [0, T] and consider the collection of random variables 
Y = col(yi, 7/2) • • • , Un) where 7/j = J Q T 4>i(t)X(t)dt. Then Y ~ iV(0, 7/v) where is the N- 
dimensional vector of all O's and In is the N x N identity matrix. 

iv) E[/ T f(s)X(s)dsX(t)] = f(t),0 < t < T and otherwise. 

This is equivalent to saying that {X(t), —00 < t < 00} is a zero mean stationary Gaussian process 
with covariance R(t, s) = R(t — s) = 6(t — s) where S(.) denotes the Dirac delta function. 
Indeed let 

R w (t)= [ S w (X)e i2wX dX = sm2nWt : -00 < t < 00 (2.3) 
J-w 

and hence for any / G L 2 we have: 

/oo 
Rw(t - s)f(s)ds = f(t) (2.4) 
-00 

where the limit is in L 2 . 

Clearly such a process is not physically realizable since by Bochner's theorem R(0) = limvi/-^oo J_ w 1-dX = 
00. in other words its sample paths cannot be in L 2 . It is worth remarking that from above: the 



3 



process Y(t) = limjy^oo Jo X\y(s)ds is a zero mean Gaussian process with variance t or is Brown- 
ian motion. The point is that the process X(t) is not well defined and thus X(t) only formally the 
derivative of Y(t). 

Herein lies the problem. Clearly we can make sense of operations when L2 functions act on X 
and in hence problems where white noise is the input to a linear time-invariant system we can give 
mathematical meaning by the limiting arguments (in q.m). However even simple nonlinear operations 
such as squaring, i.e., Y(t) = X 2 (t) do not make sense because such a process would have infinite 
mean and its covariance would be the product of delta functions that is not defined in any meaningful 
way. 

In the following section we show that the squaring operation does make sense if we perform a 
suitable renormalization of the process X\y(.) and then the limiting process itself is Gaussian white 
noise. This is indeed an unexpected result because squaring a simple Gaussian random variable results 
in a chi-squared random variable. 

3 Renormalized powers of white noise 

First note that from the definition of {Xw(t)}, —00 < t < 00 is a stationary Gaussian process, we can 
represent Xy/{t) as: 

X w(t) = §^X(0) + v w (t) (3.5) 
where vw(t) is a zero mean Gaussian r.v. independent of X(0) and variance given by: 

Let us now define the following process: 

Y w (t) = -t=(^) - = "T^ , -00 < t < 00 (3.7) 



Then Y\y(t) is a centered (mean 0) , renormalized process that denotes the nonlinear transformation 
(squaring) of the pre-white noise process. We now state and prove the main result. 



Theorem 3.1 Consider the renormalized and centered process Y\y(t) defined in 3.1 above. Then: 



lim Y w (t) = Y(t) in L 2 (P) x L 2 (3.8) 

W— ¥00 

Moreover Y{t) is Gaussian white noise. 

We prove the result through the following two results. 
Proposition 3.1 Let Ryy(t) denote the covariance ofY\y(t). Then for every /(.) £ L2 



lim 

W— >oo 



/oo 
h%r(t - s)f(s)ds = f{t), -00 < t < 00 (3.9) 
-00 



or formally: 
where 6(.) is the Dirac delta function. 



lim Rw(t) = S(t) -00 < t < 00 

W— >oo 
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Proof: 

First note that Y\y(t) is a w.s.s. process whose covariance denoted by 

E Y (t)- 1 /? 2 (t)- 1 ( sia2 * Wt Y 



Next note that by Parseval's theorem: 

I I Sill ZTTVV l \ 2 

2W 1 



1 [°° (sm2irWt\ 2 , 1 f°° o 1 



Define the measure d<pw(t) = 2W ( sm 2 Jt Wt ) 2 dt so that 4>w{— °°j 00 ) = 1 = f^° 00 4'w{'t)dt. Hence for 
any continuous /(.) G L2 we have: 



f-Jf-ooM + sWwis) - f(t)\ 2 dt = JZ\f-oc(f(t + *) - f(t))d^ w (s)\ 2 dt 

< IZI-Jf(t + *)- f{t)\ 2 dd> w {s)dt 



where we have used Minkowski's [15] inequality for integrals in the second step noting that d(f>w(-) 
defines a measure that puts mass 1 on (—00, 00). Now performing a change of variables by substituting 
2ttWs = t we can re-write: 

izf-jnt+s)- m \ 2 d^ w (s)dt = rj!°j/(t+^)-mi^(^) 2 ^ 

< -ltow(t,r)(— )||£ aXia 

where gw(t, T ) = (/(* + 2^pf) ~~ /(*)) anc ^ -^2 x ^2 = ^2(— 00,00) x L2(— 00,00) with the product 
measure defined thereon. 
Now , for every fixed r 

fZ\9w(t,r)\ 2 dt<2\\f\\ 2 
and hence using Fubini's theorem and the fact that J^^i^^-^dr = 7r 

r oo n o l^(t,r)(^)| 2 dtdr<2 7 r||/|| 2 
Moreover \gy/{t, r)| — > as W — > 00 for every r fixed we have that 

J-JJ-ooifi* + s)<hr(s)d8 - f(t)\ 2 dt -> as W -> 00 
by dominated convergence. 



The second result we prove is the convergence to a Gaussian process. For this we need the following 
result. 

Lemma 3.1 Let {Xw(t), —00 < t < 00} be a zero mean stationary Gaussian process with bandlimited 
spectral density. Let (a, 6) and (c, d) be any non- overlapping intervals in 3ft. Let h(.) £ L<i- Define the 
random variables X^ b (resp. X^) as X^ h = J a X w (s)h(s)ds (similarly for X^ d ). 
Then {X^ ,X^ d ) are asymptotically independent as W — > 00 
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Proof: To show the result it suffices to show that the random variables are asymptotically uncorrelated 
since they are jointly Gaussian by construction. 

/oo roo 
/ l^ h) (u)l {c4) (v)h(u)h(v)R w {u - v)dudv 
-oo J —oo 

Hence from applying the result of Proposition 13. II we have: 

/oo 
l (a , b) (u)l iCjd) (u)h(u)h(v)du = 
-oo 

since (a, b) and (c, d) are non-overlapping. 

Remark 3.1 From above it readily follows that the random variables YW and defined in an 
analogous way are also asymptotically independent as W — )■ oo since they are junctionals of the 
underlying process X w on non- overlapping intervals. 

Let Cw(h) denote the characteristic functional of {Yj }, — oo < t < oo defined as: 

C w (h) =TE[e l ^ W ' h H h€L 2 (0,T) (3.10) 
and [x,y] t = f£ x(s)y(s)ds for x,y G L 2 [0, T], t <T . 
Proposition 3.2 

lim C w {h) = e~^T (3.11) 

W— >oo 

or {Y^} converges to Gaussian white noise in L 2 (JP) x L 2 
Proof 

First note that from the asymptotic independence we have: 

lim n^ [YWMt+s \ = lim ne i[YW ' h]t M^ W ' h]t+S ] 

W— >oo W— >oo 

where = x(s)y(s)ds. 

Now to show that {Y t w } converges to Gaussian white noise it is sufficient to show that: 

lim dt \^ lv - = — h 2 (t) a.e. in t (3.12) 

Note that from the boundedness of e^ Y and the convergence in q.m. of Y w it is easy to show 
that: 

lim ±TE[jl YW M t] = ± Um M e ^ w ,hh] 

W^-oo dt dt W^oo 



Now: 

(3.13) 



-E[e*[ yW >]*] = lim i-E 
dt L J A^o A 



e i[Y w ,h]t + A] _ e HY w ,h]t] 



where the limit on the r.h.s. of (|3. 13[) is to be interpreted in L 2 (P). 
Hence: 

>] _ e ^ w M^ = ljE[e^ YW ^ (e^( A ) - l)] 



6 



where £ W (A) = [Y w ,h]\ +A . 

Now we make use of the following identity: 



1 + ix - ^x 2 - i I 



to obtain for each A > by asymptotic independence: 



lim JE[e l[YW ' h]t C w (A)}= lim lE[e i[YW ' h]t ]JE[(, w (A)} = 



For the second term we have: 

lim m[jP W M(Z w (A)) 2 ]= lim ne l[YW ' h]t }Vti w (A))' 

Now from Proposition 13 . 1 1 we have: 

lim E(£ W (A)) = / /i 2 (u)Ax 

For the third term we obtain: 

'0 

< (E[e M/ (A)i 2 )^ 

Therefore, once again using the result of Proposition 13.11 we obtain: 



lEfe*^'^ / s 2 e^ W ^- s Us]\ < E|^(A)| 3 
Jo 



2 



3 



t+A \ 2 



lim (E[e H/ (A)| 2 ) 5 = / h 2 {u)du) =OfA3) 
Therefore combing all the three estimates above we obtain 



e i[Y w ,h]t+A] _ gip^A]*] 



lim lim 

A-K)W/->-oo Efe ? [ r ' h l 4 



A™K(/ <+A/l(n)2 ^ + 0(Af) ) 



= h 2 (t) a.e.t 

Therefore we obtain: 

lim ne t[YWMT ] = e-^T 

or the limiting process has a characteristic functional that corresponds to standard Gauss measure on 
L2[0,T] for every T and hence is Gaussian white noise. 
This completes the proof. 



A consequence of the above result is that any integral power of white noise should remain white 
noise under proper re-normalization. Indeed it is the case and we show this below. 
Now let aw(t) = r^Io) ' anc ^ so ^ or an y integer n, using (|3.5p we obtain: 

Xfr(t) = (a(t)X w (0) + v w (t)) n = ( n )a P w (t)XZ(0)^- p (t) (3.14) 
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Prom the independence of X(0) and vw(t) we have: 

n 

RUt) = C o^(X{| / (t)^(0)] = ^a^(t)]E[X^(0)]]E[^- p (0]-(]E[X^(0)])^ 

n-l 

= a»(t)E[X#(0)] + J2^w(tMX^ P (0)M^w P (t)] - (E[X^(0)]) 2 (3.15) 

p=0 

Define: 

„ /x XZ, (t) (0)1 

YS-ft) = wl 7 1 -^tt 1 , n > 2 3.16 

where (n — 1)!! = (n — l)(n — 3)(n — 5) ... 1. 

Then we can state the following theorem about higher order powers of white noise. 

Theorem 3.2 Let {X\y(t)},—oo < t < oo be a Gaussian process whose spectral density is flat, of 
unit power and with support in [—W,W]. 

Then the process {Yfy(t)}, — oo < t < oo converges to Gaussian white noise in -^(P) x L2 as 
W — > 00. 

Proof: The proof essentially follows from the arguments in the proof of Theorem 13.11 Indeed the 
result follows from the fact that for any Gaussian N(m,a 2 ) r.v. the moment of order n is given by: 

E[(X-m) n ] = (n-l)\\<T n , (n-l)!! = (n-l)(n-3)(n-5)...l, n even 
= 0, n odd 

Note covariance of Yfy- is just given by ^ n -i^2W) n • Now from 3.15 it can be seen that R^{t) = 
C.Ryy(t) where C is a constant. 
Let /(.) G L2 and consider: 

Rn ^- s) f(s) ds = 1 r ( sm27TW{t - s) ) n f( S )ds 

Noting that < -Rw(O) we have : 

"/(-)cfa| < |/2wr(0)n|/|| < C^ll/M < 00 



(n-l)!!(2Fy)- i y_oo V 
Furthermore, see [12] : 

'sm2TrW(t - s 



dt = C^{n) < 00 



(n-l)\\{2Wy J.^ V 7r(t-a) 

Then, we can repeat the arguments as in Propositions 13.11 and 13.21 mutatis mutandis to show that 
the normalized process is Gaussian white noise whose variance depends on the constants. 



Remark 3.2 This result can be directly extended to homogeneous polynomials P(Xw(t)) in the ob- 
vious way by defining the renormalization factor to normalize the highest power of the polynomial. It 
can be shown the lower order powers do a play a role in the asymptotic limit. 
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